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Abstract 

We investigate the conjecture suggested by numerical simulations and 
experimental evidence that the scaling exponents for the Navier-Stokes 
equation are the same as for a suitable linear advection equation. We 
prove the result for the Navier-Stokes equation with additive noise in a 
2D spatial domain. To analyze the coupled system of the Navier-Stokes 
field u and the advection field w, one introduces a parameter A which 
gives a symmetric system for (u A ,io A ). This system is studied for any A 
proving the well posedness and the uniqueness of its invariant measure 
Assuming universality of the scaling exponents to the force, the scaling 
exponents of u x and w x are the same. We prove that this is true also in 
the limit as A — > 0, that is the 2D Navier-Stokes equation with additive 
noise has the same scaling exponents as the linear advection equation. 
Therefore, we have anomalous behaviour of the scaling exponents for the 
stochastic 2D Navier-Stokes equation. 



1 Introduction 

We consider the stochastic Navier Stokes equation 

du . , _ dWi 

(1) ^ ~vAu+{u-V)u = -V Pl + —, 

describing the motion of a fluid in an open bounded domain D of K d . We assume 
that T> has a smooth boundary dT> which satisfies the locally Lipschitz condition. 
Here, u — u(t, £) is the velocity vector field, p — p(t,£) is the scalar pressure 
field and v > is the viscosity. W\ — Wi(t, £) is a Gaussian random field white 
in time, subject to the restrictions imposed below on the space correlation. 
The velocity field u is subject to the incompressibility condition 

(2) v-u(t,0 = o, te[0,ruez>. 

and to the boundary condition 

(3) u(t,0 = 0, te [0,T], £ed£>; 
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an initial condition is given 



(4) u(0,0 = «o(0. CeP. 

Recently (see pQ and the references therein), it has been proposed that the 
Navier-Stokes equation and a relevant linear advection model 

, _ dw . . _ N _ 9W2 

(5) — !/Aw + (it ■ V)w = -Vp 2 



dt y L - dt ' 

have the same scaling exponents of their structure functions, even if their statis- 
tics are different (i.e., in the "jargon" of stochastic PDE's their invariant mea- 
sures are different). Here u is solution of equation (TTJ); the noises W\ and W% 
are independent and identically distributed. 

We recall the definition of structure function (see [21]): it is an ensemble 
average of power of velocity differences across a lenght scale. The longitudinal 
structure function is 

'{Ht,z+it)-u(t,o}-ir 

with £, £ + 1£ £ T>, p £ N and I £ R; £ is a versor and we take the scalar product 
of vectors in R d . Similarly one can define the transverse structure function. 

If the velocity field u is stationary in time, homogeneous and isotropic in 
space, the latter quantity depends only on the lenght scale I. We denote it by 
SP{1), that is 

(6) SS(Z) = ({[«(*, lD-u(t,0)].|}") 

for any t. For a turbulent field it is important to know how S^(l) depends on I 
for small I. The scaling exponents £ p are defined by 

(7) S*(l) cx l<* 

and are assumed to be universal in the limit v — > when / lies in the inertial 
range rj v 3 / 4 < I < l (see e.g. [2Ql [24]); that is an equivalent definition is 

is) &=i.»^. 

2^0 log/ 

Notice that I — > implies also v — > 0. 

Dimensional considerations provide Cp = §• However, numerical and exper- 
imental results give £3 = 1, £2 > § but close to the value |, and an anomalous 
behaviour Cp < f f° r P > 3. Therefore, the determination of the scaling expo- 
nents Cp is a crucial problem in turbulence theory (see [23] )■ 

However, [JJ provides numerical evidence that both the 2D Navier-Stokes 
equation and the Sabra shell model have the same scaling exponents of the 
corresponding linear advection models. If this statement were true, then it would 
allow to reduce the determination of the scaling exponents for the Navier-Stokes 
equation (fTJ) to the easier problem of determination of the scaling exponents 
for the linear advection model ([5]). As far as the linear advection problem 
is concerned, its scaling exponents show anomalous behaviour (see among the 
others, [TJ[2J[25] and the references therein). 
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It is then of interest to understand rigorously the properties of the joint 
system 



(9) 



9u . , _. _ 
— - vAu + {u ■ V)u = -Vpi -\ 

dw , . ^ 

— vAio + (it • v jiu = — Vp9 



dWi 
dW 2 



with the divergence free condition and the boundary condition. If it has a unique 
invariant measure /z, then the ensemble averages are done with respect to this 
measure, i.e. for the linear advection model ([5]) 



(10) S*(l)= ll{[w(lt[)-w(0)}-Zy(x(du,dw) 
and for the Navier-Stokes equation (|TJ) 

(11) SP(l) = U {[u(l£) - u(0)] • iV^du, dw) = I {[u(li) - u(0)] • f}"m(d«) 

being m(du) = J /i(du,dw) the unique invariant measure for ([1]). 

The analysis of system (JSJ has been performed by adding two terms: 



(12) 



du 

Th 

dw 

— — vAw + (u ■ V)w + \(w ■ \7)w 



uAu + (u ■ V)u + X(w ■ V)m = -Vpi + 



Vp 2 



dWi 
dW 2 



where A £ K is a parameter. We denote by (u x ,w x ) its solution. For A = we 
recover ((9]) and for A = 1 system (|12[) is symmetric. 

On the other hand, for any A ^ 0, system (|12p enjoys the following prop- 
erty: setting u A = Aw A and multiplying the second equation by A, we have a 
perfectly symmetric system for the pair (u A , u A ), except for the force and initial 
conditions: 



(13) 



vAu + (u ■ V)u + (v ■ V)u = -Vpi 



du 

dv 

— vAv + (u ■ V)u + (v ■ V)w 

^ at 



dWi 



-AVp 2 + A 



dt 

dW 2 



dt 



with u A (0) = uo and w A (0) = Awo- Thus, assuming the universality of the 
scaling exponents to the force, it follows that u x and v x = \w x have the same 
scaling exponents for any A ^ 0. Then the same consideration holds for the 
couple (u x ,w x ) for any A ^ 0; indeed, 5^(0 = X p S^ x (l) and © gives the 
same C, v for v x and w x . 

The crucial point is to see if this holds also in the limit as A — > 0. 

The continuous dependence of the solutions to system (fT2)l as A — > has been 
investigated rigorously in the context of certain nonlinear phenomenological 
shell model (see [1] for the deterministic case and jS] for the stochastic case). 
Moreover, [S] consider the asymptotic dynamics (for large time) on the attractor, 



lr The linear advection equation {5j has scaling exponents universal to the forcing, when 
the forces act only on large scales (this has been investigated by physicists, see e.g. |2l l9lll4| ). 
Some numerical evidence shows the same for the nonlinear equation |[TJ (see [l]). 
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proving the continuous dependence on A of the attractor. Here we address 
the problem for the stochastic 2D Navier-Stokes equation. We consider the 
stationary statistical regime described by an invariant measure. The main result 
of our paper is that the unique invariant measure /x A for system (|12p converges 
to the unique invariant measure [i for system ([9]), as A — > 0. Bearing in mind 
the expression (fTTj) . this implies that if the structure functions of u x and w x 
have the same scaling exponents for any A ^ 0, this is true also for A = 0, i.e. 
the Navier-Stokes equation ([T]) and the linear advection equation have the 
same scaling exponents of the structure functions. 

Finally, we point out that we present our results for the 2D Navier-Stokes 
equation, because of the lack of well-posedness in the 3D setting. But we do 
not use any peculiarity of the 2D case (except for the regularity results in the 
Appendix) . 

The paper is organized as follows. Section [5] gives the abstract setting and 
summarizes well known results on Navier-Stokes equation and Wiener processes. 
Classical results on the Markov semigroup are given in Section [4] and on the in- 
variant measure ^ in Section[SJ The new results are the continuous dependence 
of the solution on A (Section [3]) and the continuous dependence of the invariant 
measure on A (Section [B]). The Appendix provides regularity results, extend- 
ing the results of the previous sections from the space of finite energy to more 
regular spaces involved in the definition of the structure functions. 

2 Notations and hypothesis 
2.1 Functional setting 

We define the functional setting to study the Navier-Stokes equation. From now 
on, we consider a bidimensional spatial domain D. We refer to [TUJ [31] for the 
main results. 

Let V be the space of infinitely differentiable vector fields u on T> with 
compact support strictly contained in V, satisfying V ■ u = in V. We take the 
closure of this space in [L 2 (2?)] 2 and[_ffg (T>)] 2 ; we obtain respectively 

H = iu G [L 2 {V)] 2 ■ V -u = Q in D, u-n = on dv\ 

where n is the normal versor to dT>, and 

V = ju € [if 1 (2?)] 2 : V- u = in V, u = on <9£>j . 

They are separable Hilbert spaces with the inner products and norms inherited 
from [L 2 {V)] 2 and [H^(V)] 2 respectively: 

Hh = (w> u ) , and ( M > v ) = / [ u (0 ■ v (0] 

Jv 

\\ufy = ((u,u)), and ((u,v))= f [V«(0 • V«(0] 

Jv 

Denoting by H' and V the dual spaces, if we identify H with H', we have 
the Gelfand triple V C H C V with continuous dense injections. We denote 
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the dual pairing between u € V and v G V by (u, v)v,v ■ When v G H, we have 
(u,v) V y = (u,v). 

Let II be the orthogonal projector in \L 2 (D)] 2 onto H; then the Stokes 
operator is 

Au^-TlAu, yu e D{A) = [H 2 {V)] 2 n V. 

The operator A is a closed positive unbounded self-adjoint operator in H with 
the inverse A -1 which is a self-adjoint compact operator in H; by the classical 
spectral theorems there exists a sequence {7^}°^! of eigenvalues of the Stokes 
operator with < 71 < 72 < . . . , corresponding to the eigenvectors ej G D(A); 
{ej}j form an orthonormal basis in H. We have that jj behaves like j for 
j -> 00. 

For a > we will denote the a-power of the operator A by A a and its 
domain by D(A a ); we have |M|^ (j4q -j = Z)°1 1 7 j 2q | < u,ej > \ 2 . Moreover, 

V = DiA 1 / 2 ), and D(A ai ) is compactly embedded in D(A a2 ) for ai > a 2 . By 

interpolation we have ||m||.d(aV 4 ) — ^Mff 2 |MI v 2 ■ Finally, D(A~ a ) denotes the 
dual of D(A a ). 

Let &(•,-,■) : V x V x V — > R be the continuous trilinear form defined as 

b(u,v,z)= f ([«(£) -VM£)) 

It is well known that there exists a continuous bilinear operator B(-, •) : x 

V — ► V' such that (B(u,v), z) — b(u,v,z), for all 2 € V. By the incompress- 
ibility condition, for u,v,z G V we have (see e.g. [3T]) 

(14) (-B(w, v),z)v,v = — {B(u, z),v)v,v< and {B(u, v), v)v,v = 0. 
Furthermore, there exists a constant C such that 

(15) \\B{u,v)\\ v > < C||w|| I)(yl i/4 ) ||w|| £ , (A i/4 ) . 

This comes from Holder inequality so ||i?(tt, u)||y < ||w|| i 4||ti|| i 4, and then using 
the continuous embeddings L'(A 1 / 4 ) C [if ^(P)] 2 C [£ 4 (£>)j 2 . 

Projecting the equations of (fT2|) onto 7?, the following abstract formulation 
is obtained 




du x + [vAu x + B(u x ,u x ) + XB(w x ,u x )]dt = dWi 
dw x + [vAw x + B(u x , w x ) + XB(w x ,w x )]dt = dW 2 
u x (t ) = u Q 
w x (t ) = Wo 



with initial conditions uq, wq G H at to € R. 
The Cauchy problem is studied on [fo,oo[. 

Define H = H x H, V = V x V and L>(2 Q ) = £>(A") x D(A a ). If 1 = 
(xi, X2) G i? and y = (j/i, 1/2) G -ff, we define the scalar product in as 

{ X ,V)8 = + (^2,2/2) 

and the norms in H and V as 

Me = \ x i\h + \ x s\h, i = (n,i2)6fl 
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\\xf v = \\x 1 W v + \\x 2 \\ 2 v x = ( Xl ,x 2 )eV. 

Moreover, define the linear operator A : D(A) C H — » H, or also A : V ^ V, 
as Ace = (Aei, Ae 2 ); we have V = D(A 1/2 ). 

For every A £ K, define the bilinear continuous operator B\ from V x V to 
V as 

-Sa 0,2/) = (-B(a;i,Z/i) + XB(x 2l yi), B(xi,y 2 ) + XB(x 2 ,y 2 )j, 

where as usual we have used the notation x = (x\,x 2 ), y = (yi,y 2 ). 
By ([T ^) -([T^ ]) . we have 

Lemma 2.1 For any A 

(B x (x,y),z)yy, = -(S A (x,z),y)y (B x (x,y),y) vv , = 0. 

Moreover, there is a constant C\ > swc/i i/iai 

I|5a(^?/)||^ <Cx||x|| c(ll/4) ||y|| D( ^ /4) , Vx^e^i 1 / 4 ). 

We write (|16[) in more compact form as 

(l7] ( du x + [uAu x + B\ (u x , u x )]dt = dW, 

( ' I u x (t ) = u 

where u x — (u x ,w x ), u = (u ,wo), W = (Wi,W 2 ). 
2.2 Stochastic driving force 

The r.h.s. of the stochastic equations in (IT51) and (fTTj) are defined as follows: 

00 

(18) Wi(t)=£v©^(*) e i. 

i=i 

with (/3*-) mutually independent standard (scalar) two-sided Wiener processes 
defined on a filtered complete probability space (Q, J 7 , (.Ft), P), and qj > for all 
j. This means that W\ and W2 are two independent two-sided Wiener processes 
with covariance operator Q = diag{qj}; Q is a positive linear operator in the 
Hilbert space H. 

Let us assume that the operator Q satisfies the following property: for some 

OO 

(19) X> 7 f°<oc. 

Then Wi is a Z?(A Q °)-valued Wiener process. For details we refer to [TT]; but 
let us show the easy estimate 

3 h j 

In the following we are interested in two cases: 
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• only a finite number of qj do not vanish, i.e. 

3 n e N such that qj = V j > n. 
In this case condition (|19p is fulfilled for any ao- 

• an infinite number of qj do not vanish. In this case, the sum in (|19[) is a 

series and the sequence {qj} must satisfy qj < j- 2< *o-i-e f or some e > 
as j — ¥ oo; by this we mean 

lim qjj 2aa+1+e < oo 

(recall that jj < j and j < 7 j as j — > 00). Notice that if condition (|T9^) is 
fulfilled for a value «o then it holds for any value smaller than ao- 

In the sequel, we consider ao > — \ in (fH?l) . If ao = — j then the processes 
W x and W 2 have paths in C(K; D(A^ 1 ^ 4 )) P-a.s. If a = then Q is a 
trace class operator, i.e. TrQ :— Ylj=i Qj < 00 • 

The case of finite dimensional noise is the interesting one in turbulence theory 
(see [24]) and in numerical simulations: the force is injected at large scales (j 
small) and is transported by the non linearity to small scales (j large). However, 
the case of non degerate (or full) noise, i.e. qj 7^ for all j, is interesting from the 
mathematical point of view, in the sense that for some hydrodynamical models 
the unique invariant measure for system (|13D is equivalent to the Gaussian 
invariant measure of the associated linear Ornstein-Uhlenbeck system. This has 
been proved in [18 for hyperviscous fluids, and can be proved for shell models 
following the setting of [31 [12] . This may be useful to get concrete expression for 
the invariant measure for system (|13p (dealing with hyperviscous fluids or shell 
models). By the way this could help in checking the universality of the scaling 
exponents to the force. However, the results of [18] give so far the expression 
of the Radon-Nykodim derivative, obtained by Girsanov transform, only for the 
solution process on any finite time interval and not for the stationary process 
(even if a formal expression is given in [25) for Navier-Stokes and applies also 
to our model (flU)0 . 

3 Continuous dependence of the solution on A 

We begin giving the usual definition of solution of (|17p . This is a strong solution 
in the probabilistic sense and a weak (or generalized) solution in the sense of 
PDE's. 

Definition 3.1 We are given a stochastic basis (f2, J 7 , (J-j) , P), a D(A^ 1 ^ 4 )- 
valued Wiener process W , a time interval [iojT] and o- n initial data uq £ H. 
We say that a stochastic process u x is a solution of the Cauchy problem j 17\ ) on 
[to,T] if it is a continuous Tt-adapted process in H, 

u x e C([t , T]; H) fl L 2 (t , T; Z^l 1 / 4 )) P - a.s. 
and (|17[) is satisfied in the following sense: 

(20) (u x {t), V) + f iy{A^ 4 u x {s), A 3 / 4 ip)ds - [ (B x (u x {s), i>) , u x {s))ds 

J t J to 

= {u , V) + (A^^Wit), l^V) - (A- 1/4 W{t ), A 1 / 4 ?/*) P - a.s. 
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for any t G [t ,T] and i/j G D{A 3 / 4 ). 

Recalling Lemma \2. II we notice that 

(B x (u\u x ) » - -(B X (u\^j),u x ) 

and the latter expression is well defined. Therefore equation (|20p corresponds 
to (g7j). 

The well posedness of the stochastic Navier-Stokes equation is a well estab- 
lished result; we cannot give a complete list of references (see e.g. references in 
|20j). But we refer to the results proved for the 2D case with additive noise. 
Indeed, following [2TJ[T7], one can prove existence and uniqueness of the solution 
to system (fT7)l . The procedure is the same as for the Navier-Stokes equation, 
since the operator B\ enjoys the same properties of the bilinear operator B of 
Navier-Stokes, as given in Lemma \2. II 

Theorem 3.2 Let A G R. Let us assume that uq € H and that assumption 
(|19p is satisfied with «o = ~\- Then, there is a process u x with paths in 
C([to,T]; H) n L 4 (to, T; D{A l / A )) P-a.s., which is the unique solution u x for 
(|17p in the sense of Definition lS.R 

In the proof, one introduces the Ornstein-Uhlenbeck process z solution of the 
(linear) Stokes type problem dz + vAz dt = dW, z(to) = 0: 

z(t)= f e-" A(t - s) dW{s). 

We have (see [TT] ) 

Proposition 3.3 Let assumption (|19[) be satisfied. Then, the process z is a well 
defined Gaussian, ergodic and continuous process in D(A a ° + ^). In particular, 
for any finite interval \t$,T\ there is a P-a.s. finite random variable C = C(lo) 
such that 

(21) sup ||2(t,w)|| j ao+i <C(u) forP-a.e.Lu. 

t <t<T Ul - A > 

In particular, for ao = — \ we have z £ C(R; D^A 1 ^)) P-a.s. 

The aim of this section is to prove the continuous dependence of the solution 
u x on A. 

Theorem 3.4 Let assumption (|19p be fulfilled with ao = —4. Then, for any 
Ao G R, Mo € H and finite time interval [to,T], we have 

lim sup \u x (t) - u Xa {t)\~ = P-a.s. 

A ^ A o t <t<T 

Proof. We prove the theorem only in the case Ao = and £o = 0, the general 
case being the same. Since ao = — \ we have that z G C([0, T]; D(A 1 / 4 )) P-a.s., 
and there exists a random variable C (which is finite P-a.s.) such that 

(22) ll z llc([o,T] ; _D(I 1 /4)) < C P-a.s. 
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Step 1 (preparation). Let us define the new process 

q x := u x + Xw x 
and the corresponding difference 

p X :=q X -q =q X - u a . 

Since u x and w x satisfy 

du x (t) + [vAu x {t) + B (u x (t),u x (t)) + XB (w x (t) , u x {£))} dt = dWi{t) 

dw x (t) + [vAw x (t) + B (u x (t) ,w x (t)) + XB (w x (s) , w x (t))] dt = dW 2 (t) 
then q x and p x are solutions, respectively, of 

dq X (t) + [vAq X {t) + B (q X (t),q X (t))} dt = dWi{t) + XdW 2 {t) 

dp X {t) + [vAp X (t)+B (q X (t),p X (t))+B (p X (t),q X (t))+B (p X (t),p X (t))} dt = XdW 2 (t). 

These equations are obtained by using the bilinearity of the operator B and 
have to be understood in the integral sense (in time) and with scalar product 
with suitable text functions, as in Definition 13. f I 

Step 2 (bound on q x ). Let us prove next that P-a.s. 

(23) sup sup |g A (t)| < oo, 

AG[-1.1] 0<t<T 

and 

(24) sup f \\q X {t)\\ A D(A1/4) dt<^. 
Ae[-x,i]Jo 

Notice that the equation satisfied by q x is the Navier-Stokes equation with initial 
data uq + Xwq and forcing term dW\ + XdW 2 ■ Therefore we get a priori estimates 
following the technique of [21]. We define 

v X := q X - z X , 

where z x is the solution of the linear problem 

dz x + vAz x dt = dWi + XdW 2 , z A (0) = 0. 

In this way we get rid of the stochastic force. Indeed v x satisfies P-a.s. a 
deterministic equation (no stochastic integral appears): 

^V(i) + vAv x {t) + B (v x (t) + z x (t),v x (t) + z x (t)) = 0. 

As we shall see, the paths of the process v x are more regular than those of q x 
and z x . We are going to prove that v x e C([0, T];H)C\ i 2 (0, T; V) P-a.s.. From 
now on in the proof we proceed pathwise. 

Formally, taking the scalar product in H with v x and using (jT4")) we get 

1 f 1 

^\v X {t)\ 2 H + v \\v x { S )\\ 2 v ds <-\u + Xw \ 2 H 

+ f \(B(v x (s) + z x (s),v x ( S )),z x ( S ))\d S . 
Jo 
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Rigorously, the above inequality can be proved either by general abstract theo- 
rems or by taking finite-dimensional approximations, performing the computa- 
tions at the finite dimensional level and then taking the limit (see |31j). Then 
we have 

(25) 

1 I \, . i2 f n \ . N ||2 1 . ,2 

-\v {t)\ H + v \\v {s)\\ v ds- -\u q + \wq\ h 
< f I (B(v x (s) + z x (s), « A (s)), z x {s)) I ds 



<cy o |u A (*)i^ /2 i|i; A (*)ii^ a ii« A wib ( xv*)d*+c y o ii^( s )ik||z A ( s )||; (A1/4) 



ds 



< 



From now on, C„ denotes a constant depending on v 1 which may vary from line 
to line. Hence, 

tA^\|2 ^ i, _, \„,, _|2 i IUA||4 / I-A/M2 . , ^ ,|i A||4 



|« < K + Awo|# + C*, y ||z || C([ o, T ] ;r ,(Al/4)) / \ v {s)\ H ds+C u t\\z || C ( [0 ,T];£>(AV4))- 

J 

By Gronwall lemma we deduce 

(26) \v x {t)\ 2 H < (\u + \w \ 2 H + l) e C " l|2A|| c([o,T, iDCA i/4„' Vte [ 0)T ]. 
thus 

sup sup |v A (i)|# < oo P — a.s. 

Ae[-l,l] 0<t<T 

Since g A = v x + z A , this implies ([2"3"f . using the bound (f2"2"|) on the paths of z x . 
From ([25]) and <[26j) we get 

sup / \\v x (t)\\ydt < oo, 
Ae[-i,i] Jo 

and therefore 

sup / ||¥ A (t)||4 1/4) cft< ( sup sup |U*(i)|?r J ( SU P / ll« A (*)llv*] <oo. 

AG[-1,1]J0 \AG[-1.1] 0<t<T J yAG[-l,l]JO / 

From this ([MJ) follows. 

Step 3 (convergence of p x ). Next we prove that 

(27) lim sup |p A (<)| =0, 

A ^°0<t<T 

and 

(28) lim^ T ||p A (t)||^ 1/4) dt = 0. 

In one sentence, this is a consequence of the various bounds that we have 
established previously, and the fact that the initial condition \w$ and the forcing 
term \W2if) converge to zero, as A —> 0. 
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As before, define the difference 

v x := p x - Xz, 

where 

z(t)= f e- vA{t - s) dW 2 {s) 



o 



is the solution of dz(t) + vAz(t)dt = dW2(t), z(0) = 0. We work on v x in a 
similar way we did for v . Then 



d X (t) + uAv X (t) + B (q x {t), p x (t)) + B (p x (t), q x (t)) - B (p x (t), p x (t)) = 0. 



d 
~dt 

By (|f4l) and (|T5|) . taking the scalar product with v x in H we have 
\{B (q X ,p X ),v X )\ = \(B {q X ,v X ),p X )\ = \(B (q X ,d X ),Xz)\ 

< c||9 a IId(ai/4)||v a ||v II a^|1 jD(A i/4 ) 



< g ll^llv + aAl|g A ||^ (AV4) + C u \'\\z\r D(A i,<y 
In a similar way we estimate 

\(B(p X ,q X ),d X )\ = \(B(p X ,d X ),q X )\ = \(B(v X + Xz,d X ),q X )\ 

< \(B(d X ,d X ),q X )\ + \(B(Xz,d X ),q X )\ 

< g ll wA llv + c^lk A lli)(A i /4)l wA ll/ + c l yX 2 \\q x \\ 4 D ( A1/i) + c v x 2 \\z\\ A D ^ Al/i ^ 

and 

\(B(p x ,p x ),d x )\ = \(B(p x ,d x ),p x )\ = \(B(v x +Xz,v x ),Xz)\ 
< \(B(v x ,v x ),Xz)\+X 2 \(B(z,v x ),z)\ 



6 

Hence 



< ^ll« A ||y + CA 4 \\z\\ 4 DiA ^ } \v X \ 2 H + C^\\z\\% (A1/i) 



\v\t)\ 2 H + v j \\v\ S )\\ 2 v ds < A 2 KI^ + ay o (11^11^/4) + A 4 ||z||2, (A1/4) ) \v X ( S )\ 2 H d S 

+ C^A 4 ||z||^ A l/4- ) + CyX 2 {^(^Wd^A 1 / 1 ) + H Z llz5(yl 1 /'4)) • 

As usual, by Gronwall lemma we get 



lim sup \v x {t)\ H = 0; 

A - >0 0<«<T 



in addition 

lim / ||u A (s)||^ds = 0. 

As in the previous step, they imply the claims ([27)1 and (|28p . 

Step 4 (convergence of w x ). We define £ A := w A - w°, and prove that 

(29) lim sup \£. x (t)\ H = 0, P - a.s.; 

o<t<T 
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this is the statement of the theorem for the second component of u . 
We have 

j t Z x {t) + vA£, x {t) = B(u°(t),w°(t)) ~ B(u x (t),w x (t)) -XB (w x (t) , w x (*)) 
= B (q (t),w (t)) - B (q x (t),w x (t)) 
= -B{p x (t),w°(t))-B( q x (t),Z x (t)). 
As in the previous steps, we deduce that 

\e(t)\ 2 H +v tweeds 

JQ 

< rKB(/( S ),^( S )) )U ; ( S ))|d S 
Jo 

<\\ WemUs + CvJ^ ||/(.s)||^ 1/4) ||u; ( S )||^ (A1/4) d S . 

Hence, 

Using (gSJ we get (ggj) . 

Step 5 (convergence of u x ). We simply notice that 

\u x (t) - u°(t)\ H = \q x (t) - q°(t) - Xw x (t)\ H = \p x (t) - Xw x (t)\ H ; 

therefore, by the results of steps 3 and 4, we have 

lim sup \u x (t) - u°(t)\ H = 0. 

0<t<T 

The proof is complete. □ 

4 Transition semigroup 

Let us denote by u x (t; x) the solution of (fTTj) with initial condition x, by Bb{H) 
the space of Borel bounded functions <p : H — > K and by (i? ) its subspace of 
continuous bounded functions. Define the operators P x : Bb{H) — > Bb{H) as 

(30) (P X ip)(x)=E[<p(u x (t;x))}. 

The process u x is said to be a Markov process in H if 

(31) E[^{u x {t + s;x))\F t ] = (P x <p)(u x (t;x)) P - a.s. 

for any y> G Cb(H), t,s > 0, x E H. Taking expectation in (|3ip one obtains the 
semigroup property P x +S = P X P X on Cb(H). 

The Markov semigroup {P x }t is said to be Feller in H if 

F f A : C b (H) -► C 6 (ff) Vt. 
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There is an extensive literature for the 2D Navier-Stokes equation driven by an 
additive noise. In particular we refer to 2U], Section 3.4, for the proof of the 
Markov property and to [21] for the proof of the Feller property. Notice that 
the proofs of [201 EH] for the 2D stochastic Navier-Stokes equation apply also to 
our system (|17|1 . Therefore we have 

Proposition 4.1 With the same assumption of Theorem \ 3.2l the solution u x 
of system (|17p defines a Markov semigroup which is Feller in H . 

As a consequence of the continuous dependence of the solution u x on A given 
in Theorem 13.41 we have the following result. 

Corollary 4.2 With the same assumptions of Theorem \3-4\ for any £, Ao G M 

we have 

(32) lim P t VO) = P?°<p(x), y<peC b (H), xe H. 

A-s-Ao 

Proof. Keeping in mind definition (|30p . this is a consequence of Theorem 13.41 
and the Lebesgue dominated convergence theorem. □ 



5 Invariant measure 

A probability measure /i A is invariant for P x in H if 

(33) fp t x <p(x)dn x (x)=f<fi(x)d^ x (x), Vtp£t3 b {H), Vt. 
Jh Jh 

We restrict to (f € Cb(H), since the Markov semigroup is Feller. 

For any A £E K one can prove existence and uniqueness of the invariant 
measure for system (fl7)) . following the lines of the proofs for the 2D Navier- 
Stokes equation. We briefly recall the results. 

It is not difficult to prove the existence of an invariant measure for the 
semigroup P x . Indeed, existence is proved assuming again (1191) with = — j 
(see [21] Theorem 3.1). For the sake of simplicity we sketch the proof of the 
existence following [8., which requires the operator Q to be trace class, i.e. 

CtQ = 0. 

First, use Ito formula on the process |w A (i)|^ 

d|w A (i)U = -2v\\u x (t)\\}dt- 2(B(u x (t),u x (t)),u x (t))dt 

(34) L. 

+ 2 <u x (t),dW > +TrQ dt, 

where TrQ = 2TrQ. Then use (fT4| . integrate in time and take expectation 

(35) IE|« A (i)|| + 2v [ E\\u x (s)\\ 2 ~ds = \x\% + 2(TrQ)t. 

Jo 

Define the family of measures {to^}t>o on Borel subsets r of H 
m^{T) ~yJ q P{u X (t;0)eT}dt. 
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Since (|55|) implies 

2v [ E\\u x {t;0)\\ 2 ~dt < 2{TrQ) T , 
Jo 



then 



(36) J \\x\\ 2 ~dm^(x) < foranyT>0. 

By Chebyshev inequality 

1 TtCI 

m x T {x : \\x\\l >R)< 



R v 

that is the family of measures {iti^}t>o is tight in any space which contains 
V as compact subset. Let us say that it is tight in H. Now, thanks to Feller 
property we can apply Krylov-Bogoliubov method to get that there exists at 
least one invariant measure /z A for system (jTTJ) (see e.g. [121 HO])- This /x A is the 
weak limit of a subsequence of ; moreover l|36p implies 

(37) / \\x\\U^(x) < T ' Q 



v 

Remark 5.1 Inequality (|37p is indeed an equality, i.e. 
(38) / \\x\\U^(x) TrQ 



v 



for any invariant measure fi obtained as the weak limit of a subsequence of 
{m^} in the Krylov-Bogoliubov method. Indeed, let u* t be the stationary process 
solving (|17[) and whose law at any fixed time is fi x . Then (|35p becomes 



E\u* t (t)\s + ^ / n^MVyds = E\u^(0)\~ + 2(Tr Q)t 
Jo 

and p7[) implies that each term is well defined. By stationarity, we get 

2v f E||^ t (s)|||ds = 2(Tr Q)t, Vi > 0; 
Jo 



hence 

The question of uniqueness of invariant measures is more subtle. In the 
last twenty years many results on uniqueness appeared; see, among the others, 

H Hao E31 123 123 EDI- 

The best result of uniqueness of the invariant measure has been proved by 
means of the irreducibility and asymptotic strong Feller property by Hairer and 
Mattingly |27) assuming that the noise acts on first few modes, i.e. 

^ H Q f ^ 3m G N, 3ji, j 2 , ...,J m eN such that 

Iji ) 9j2! • • • ) Qjm 7^ an d the remaining qj =0. 

Actually the number m and the indices j\,j2, ■ ■ ■ ,jm have to be chosen in a 
proper way, but the choice is independent of the viscosity. In particular, one 
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can take m = 4 when the spatial domain I? is a box and periodic boundary 
conditions are assumed (see general domains can be considered in the 

same way (see [25]). 

On the other hand, the first result of uniqueness of invariant measure for 
the stochastic 2D Navier-Stokes equation has been proved by Flandoli and 
Maslowski [23 assuming irreducibility and strong Feller property (see |12j for 
this technique for general SPDE's). This requires the noise to be full (qj ^ 
for all j) and with suitable space regularity 

3e > : D{A 1 /' 2 ) C R(y/Q) C D(Ai +e ) 

R(y/Q) being the range of the operator y/Q. Since D(A 1 ^ 2 ) C R(^/Q) is in- 
compatible with TrQ < oo, we write the assumption on the covariance operator 
Q in the more general form given in 15 , which allows Q to be a trace class 
operator: 

(HQi) 3e > andpe]i,i[ : D(A 2p ) C R(y/Q) C D{Ai+% +£ ). 

For instance, one can take qj = Cj~ a with 1 < a < 2. 

From now on, we assume QHQf| ) or ( |HQi[ ) in order to get the uniqueness of 
the invariant measure fi X for system (TIT)) . 

Combining the existence and uniqueness results we get 

Theorem 5.2 Assume Tr Q < oo. 

If ( HQf ) or ( |HQi[ ) hold, then for any A £ K system fj 1 T[) /ias a unique invariant 
measure Moreover it is ergodic, i.e. 

(39) lim - / P t x ip(x)dt = [ <p d/x\ V<p G C&f-ff) and Vx G 
and 

V V 

Finally, we have that the family {/^ A },\£B of these (unique) invariant mea- 
sures is tight in H , that is 

(40) Ve > 3K e compact subset of H : inf fi x (K e ) > 1 - e. 



Theorem 5.3 With the assumptions of Theorem 1 5. ,21 we have that the family 

Proof. Let us consider (|3"5)) and denote by Kr = ja; € V : ||a;||^<i?| and 
by its complementary; J^r is a compact ball in H . Then 



IHI?- 



Tr Q 1 



R n v 7 ~ f i? 
Tightness follows choosing R — in (|40l) . □ 
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6 Continuous dependence of the invariant mea- 
sure on A 



In this section we prove our main result on the invariant measures \i 



\ 



Theorem 6.1 With the assumptions of Theorem ] 5.0$ we have that /i A is weakly 
convergent to /i A ° , whenever X converges to Xq . 

Proof. We prove the theorem only in the case Ao = 0, the general case being 
the same. We previously proved that the family {^ a }a is tight in H. By 
Prokhorov's theorem, any sequence {^ A<M }fceN with lim^oo A( fc ) = has a 
subsequence {fi x 3 }jeN weakly convergent to some measure j3: 

(41) lim ftp d^"^ = ftp d(3 V<peC b (H). 

i-*x> J J 

We need to prove that j3 = /jP independently of the chosen sequence. 

By using the fact that the measure /i A and /i° are the unique invariant 
measures with respect to P A and P ( ° respectively and (1591) holds, then for any 
(p G Cb(H) we have for any x G H 

i r T 

ip dfj, X = Jim - / P t x V (x) dt 



H T ^ 



lim ± f (P A ^(x) - P<V(*)) dt + lim i f P?<p(x) dt 



= £ A + f ip d(x°. 
Jh 



Now, we use Corollary |4.2l and that for any A G R and x G H we have \Pf<p(x)\ < 
||y||c b - By means of the dominated convergence theorem, we infer that for any 
x G H 

i r T 

lim £ A = lim lim - / (P t X p(x) - P?<p(x)) dt 

A^O A^O T^oo 1 Jq 

i r T 

= lim lim - / (P t V(^) - P?<P(?j) dt 

T— >oo A— *-0 1 Jq 
rT 



= lim i / lim {P^{x) - P t V(x)) dt 

T->OC 1 Jq A^O 



= 0. 

Now, recalling (fill we get that 

(42) / Lpdj3= I ipdfi 

Jh Jh 

for any arbitrary <p G Cb{H). Thus, (3 — □ 
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A Regularity results 

We have proved our results working in the space H of finite energy. But one 
needs the velocity vector field u to be regular enough in space and p-integrable in 
order to define the structure functions (jTUJ) - (JTTJ) . Here we summarize regularity 
results, providing that each measure /i x has support contained in [C(2?)] 2 and 
every p-moment is finite. 

In the previous sections the results on the continuous dependence of the 
solution and of the invariant measure on A have been proved in H, but can be 
extended in the same way in more regular spaces. 

Following Appendix 1 of [53] we have 



Theorem A.l Let uq be a J-o -measurable random variable and p > 2. 
i) //E[|u ||] < co an 
system (|17[) such that 



i) 7/E[|uo|^] < oo and TrQ < oo, then there exists a unique solution u x for 



E 



sup \u\t)f 

0<t<T " 



"i'E [ T |S A (i)|^ 2 ||£ A W||^ < C(p,T 7 E\u\0)\l,TrQ). 



ii) Let the spatial domain T> be a torus. 

7/E[||uo||~] < oo and Tr(AQ) < oo, then there exists a unique solution u x for 
system (|17p such that 



E 



sup \\u x (t)\\ P y 

0<t<T v 



-vpE I \\u\t)\\t- 2 \\u\t)\\l {A) dt < C(p,T,E\\u\0)\\t,Tr(AQ)). 



o 



The basic tool to prove this result is Ito formula on the process \u (t)\ 



A, A IP 



<«|S A (*)|^ — ^pi^c*)!^ 3 !!^^)!!^ - p|« A (*)|^- 2 <s A (« A (*),« A (*)),« A (*)>£ft 

+ P\u\t)\ P n 2 < u\t),dW > +\p{p-l)\u\t)\v~ 2 2TrQ dt. 
or on the process ||u A (£)||^ 

+ P\\u X (t)\\ P y 2 < Au x (t),dW > + l -p{p~l)\\u\t)f~ 2 2Tr{AQ) dt. 

The latter relashionship is studied by using (B\(x, x),Ax) = which holds 
only when the spatial domain I? is a box and periodic boundary conditions are 
assumed; indeed, this property comes from the conservation of vorticity for an 
inviscid fluid. 

A consequence of this theorem is that for any deterministic initial data in 
H , when TrQ < oo for the solution u x one has 



sup \u\t)\ p 

0<t<T . 



E 

Therefore for any invariant measure 



n ,E j T \u\t)\ p ~ 2 \\u\t)\\l 7 dt<oo. 







\x\ p B d[i x (dx) < oo, / \x\ p ~<\\x\\ p ~r z d^{dx) < oo. 

H I H V 
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In a similar way we work under the assupmtion Tr(AQ) < oo. 

Finally, we generalize the result for the invariant measures. 

Theorem A. 2 Let p>2. 

i) IfTrQ < oo, then 

J \x\^\\xf^(dx) = ^(TrQ) J |z||-V(^). 

ii) IfDisa torus and Tr(AQ) < oo, then 

J M^ 2 \\4l (I f\dx) = P-^(Tr(AQ)) J \\x\ff 2 ^(dx). 

To prove this result it is enough to work on the stationary process whose law at 
fixed time is /i A and use the previous identities from ltd formula. 

We point out that when the operator Q is full, the uniqueness of the invariant 
measure can be proved under the assumption Tr(AQ) < oo using the results of 
[To] (a further generalization of ( |HQi[ )). 

Item ii) with p = 2 of the latter theorem implies that 

/ \\x\\ 2 D[X) n\dx) < oo; 

so the measure /i A has support in D(A). 

We conclude working on the torus, to get the p-integrability of the [CiV)] 2 - 
norm. When I? is a torus and Tr(AQ) < oo, let C 7 be the constant such that 
\\x\\y < Cy\\3>\\ jjr£y Then we have 

/ M P ~ 2 \\x\\ 2 D{lf \dx) < { l-^Tr{Q)[Tr{AQ)]^C^ 

for any p even and greater than 2. This is proved by induction using the re- 
lashionships of Theorem IA.2I Now, by interpolation, for any p > 2 we have 
IM| P + 1 < C v \\x\T 2 ~i then usin § that H 1+ \V) C C{V) for any 

e > 0, we get 

/ \\x\\\ c(v]]2 p\dx)<C p J ||x||f 2 ||a:||^ ) M A (da ; )<oo. 
Therefore the integrals appearing in the structure functions 

^(0 = / {Mb ~ Mo)] ■ iV^idx), 

s p w ,(i) = J{[x 2 (ii)-x 2 {o)]-ir fJ ,\dx) 

are well defined for any p > 2 when the spatial domain I? is a torus and 
Tr(AQ) < oo. 
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